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We investigate the long-range phase coherence of homogeneous and trapped Bose gases as a 
function of the geometry of the trap, the temperature, and the mean-field interactions in the weakly 
interacting limit. We explicitly take into account the (quasi) condensate depletion due to quantum 
and thermal fluctuations, i.e., we include the effects of both phase and density fluctuations. In 
particular, we determine the phase diagram of the gas by calculating the off-diagonal one-particle 
density matrix and discuss the various crossovers that occur in this phase diagram and the feasibility 
of their experimental observation in trapped Bose gases. 
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I. INTRODUCTION 

The weakly-interacting Bose gas in three dimen- 
sions has been studied in great detail over the past 
50 years [jj. Below a critical temperature T c ~ 
(2Trh 2 /mk B ) [n/C(3/2)] 2/3 , where n is the total density, 
the gas is in a Bose-Einstein condensed state. As with 
almost all thermodynamic phase transitions in three di- 
mensions, a Bose-Einstein condensate leads to long-range 
order in the system. In this particular case, the long- 
range order is determined by the behavior of the off- 
diagonal one-particle density matrix (ip'' (x)ip(0)}. If the 
one-particle density matrix goes to a constant n c in the 
limit |x| — ► oo, a Bose-Einstein condensate is present and 
n c is the condensate density. 

The physics of one and two-dimensional Bose gases 
is very different from that of the three-dimensional one, 
which makes these low-dimensional systems very inter- 
esting. From a theoretical point of view, this difference 
is caused by the enhanced importance of the phase fluc- 
tuations (3-EJEH3 In fact, the phase fluctuations are so 
large that in a homogeneous one-dimensional Bose gas at 
all temperatures and in a homogeneous two-dimensional 
Bose gas at any nonzero temperature, Bose-Einstein con- 
densation cannot take place. This is the content of the 
Mermin- Wagner-Hohenberg theorem |y, Q • 

A natural question then arises: What happens if we 
have a large three-dimensional box with a Bose-Einstein 
condensed gas and squeeze two of the sides so as to ob- 
tain a one-dimensional system? Similarly, what happens 
if we take the same box containing a Bose-Einstein con- 
densed gas at a nonzero temperature and squeeze one 
of the sides so that a two-dimensional system results? 
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These are questions that we would like to address in 
the present paper. Until recently, these questions could 
not be discussed on the basis of a microscopic theory, 
since no accurate equation of state existed for one and 
two-dimensional Bose gases. However, we have recently 
developed an improved many-body T-matrix theory for 
partially Bose-Einstein condensed atomic gases by treat- 
ing the phase fluctuations exactly H . This mean- field 
theory is valid in arbitrary dimensions and accounts also 
for the (quasi)condensate depletion. It is therefore very 
suitable for describing the crossovers mentioned above. 

Low-dimensional Bose gases are presently receiving a 
large amount of attention due to the recent experimen- 
tal realization of one- and two-dimensional condensates 
in traps llCtllll as well as one-dimensional gases on mi- 
crochips |l2t Il3l IT^ | . Because these experiments deal 
with trapped, and therefore inhomogeneous, Bose gases 
it is important to consider this situation as well. The 
Mermin- Wagner-Hohenberg theorem 0, Q theorem is 
valid only in the thermodynamic limit, and so does not 
immediately apply to trapped Bose gases. In the trapped 
case, it turns out that under certain conditions the phase 
is coherent over a distance of the order of the size of the 
system and a "true" condensate is present. Under other 
conditions, the phase is coherent over a distance less than 
the size of the system and only a so-called "quasiconden- 
sate" G3 is present @ H, H Q. In fact, this situation 
can even occur in e long ated three-dimensional Bose gases 
El Elll E3 HI ET As a result, similar questions as 
above arise: What happens to the phase fluctuations in a 
one-dimensional trapped Bose gas as we vary the temper- 
ature and the trap frequency? This is another question 
that we consider here. 

The paper is organized as follows. In Sec. [HI we 
briefly discuss the modified Popov theory presented in 
Refs. ii. In Sec. IIIII we discuss dimensional crossovers 
and finite-size effects in a homogeneous Bose gas. In 
Sec. IIVI the temperature crossover in both homogeneous 
and trapped Bose gases in one dimension is examined. 
Finally, we summarize and conclude in Sec. IVl 
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FIG. 1: Illustration of the geometry for a) the dimensional 
crossover from a three-dimensional to a one-dimensional 
Bose gas, and b) the dimensional crossover from a three- 
dimensional to a two-dimensional Bose gas. 



II. MODIFIED POPOV THEORY 

It is well known that the usual Popov theory for par- 
tially Bose-Einstein condensed gases suffers from infrared 
divergences in the equation of state at all temperatures 
in one dimension and at any nonzero temperature in two 
dimensions. These infrared problems arise because only 
quadratic fluctuations around the mean field have been 
taken into account. Specifically, the annihilation opera- 
tor for the field is written as i[>(x) = ^/uq + ^>'(x) and 
terms in the Hamiltonian that are of third order or higher 
in the fluctuating field $'(x) are neglected. In Refs. 0,0, 
it was shown that the problem of infrared divergences in 
the equation of state can be solved by taking into account 
phase fluctuations exactly and not only up to second or- 
der. The quadratic contribution to the density from the 
phase fluctuations is no(x( x )x( x )), while an exact result 
would give no contribution at all. This follows from the 
identity ^(e'^We^W) = n (1 + (x(x)x(x)) + ...) = 
1. In order to obtain the correct result, we must there- 
fore subtract the quadratic contribution from the phase 
fluctuations. This ultimately results in the following ex- 
pressions for the density n and chemical potential u 
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hu; k n T 2B (-2/i) 
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(1) 



(2n - 7i )T 2B (-2u) = (2n' + n )T 2B (-2u) , (2) 



where no is the (quasi)condensate density and n' = n—no 
represents the depletion of the (quasi)condensate due to 
quantum and thermal fluctuations. It is important to 
point out that this theory takes the (quasi) condensate 
depletion into account, and can thus work even in the 
limit ofno/n <C 1, as will become apparent in subsequent 



sections. The Bogoliubov quasiparticle dispersion rela- 

1/2 

tion is given by tiui^ = [e 2 . + 2n T 2B (— 2u)ek] , where 
£k = fi. 2 k 2 /2m. Note also that the energy argument of 
the T-matrix is —2fj,, because this is precisely the energy 
it costs to excite two atoms from the condensate |22l l23| . 

The Bose distribution function is given by N(e) — 
l/[exp(e/fcsT) — 1], where ks is the Boltzmann constant 
and T is the temperature. In addition, T 2B is the two- 
body T-matrix and V is the volume of the system. The 
two-body T matrix takes into account successive two- 
particle scattering in the vacuum, but neglects many- 
body effects of the medium. These many-body effects 
can be taken into account by using the many-body T 
matrix instead of the two-body T-matrix 9]. For sim- 
plicity we shall use the two-body T-matrix for the most 
part of this work. However, in Sec. IV B we will actually 
use the many-body T-matrices to make the calculations 
for the trapped case as accurate as possible. 

It is important to realize that the momentum sum on 
the right-hand side of Eq. ||TJ includes the term with 
k = 0. This is a result of the fact that if we include fluctu- 
ations we have that (tpQtpo), where tpo = J dx(^(x)/V), 
is not exactly equal to no. This difference is most 
important at nonzero temperatures and by taking the 
limit k ► in Eq. ft), we find that the contribution 
from the zero-momentum state to the density equals 
n + (l/V)(fc B T/2n T 2£! ). This result can be physically 
understood by realizing that (tpQ^o) can be found by cal- 
culating the average of \ipa\ 2 with a probability distribu- 
tion that obeys 



P(M) = e~ 



-(-n |i>o| 2 + 3|i/'o 



(3) 



To be consistent, the above calculation has to be per- 
formed in the Bogoliubov approximation as well. 

We define the normalized first-order correlation func- 
tion <?(x) by 



(^(x)^(O)) = v/n (x)n (0) g(x) . (4) 
In the large- |x| limit, it thus takes the form 

g( X ) ~ e -^[xW-x(0)] 2 ) _ (5) 

For a homogeneous Bose gas, the phase fluctuations are 
determined by H, |jj 
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[1 — cos(k-x)] 



(6) 



The analogous expression for the one-dimensional har- 
monically confined condensate has also been calculated. 
In the Thomas-Fermi limit and neglecting quantum fluc- 
tuations, it reads 



ix(z)-x(o)Y 



AttkI 
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B 



1 - (z/Rtf) 2 



(7) 



where Aj = y/(J + l/2)/x'/?k^, B 3 = ^C? + 1/2)7^//*', 
and Pj(z) is the Legendre polynomial of order j. Here 
fi' = fi — 2kti'(0), ra'(0) is depletion in the center of 
the trap, and the coupling constant k — a/2ir£ 2 L re- 
sults from averaging the three-dimensional two-body T- 
matrix, which is equal to Airaft 2 /m with a the s-wave 
scattering length, over the cross-sectional area of the 
one-dimensional condensate. Furthermore, £± is the har- 
monic oscillator length in the radial direction. In the ax- 
ial direction, the trap length is l z = y/K/ mu> z where w z 
is the harmonic oscillator frequency in the z-direction. 
The quantity Rtf denotes the temperature-dependent 
Thomas- Fermi radius, defined by the point at which the 
quasicondensate density n (z) = (//' — V tlap (z)) van- 
ishes. Here V tr&p (z) is the harmonic trapping potential 
in the axial direction. The frequencies u>j are given by 

= VJti + !)/ 2 

If the function <?(x) aproaches a constant in the large- 
|x| limit, the system contains a true Bose-Einstein con- 
densate. If g(x) goes to zero algebraically in the limit 
|x| — > oo, the system is said to contain a quasiconden- 
sate. Finally, if c/(x) goes to zero exponentially fast, the 
system is in the normal state. In the following, we calcu- 
late the correlation function g(x) for various situations. 



III. DIMENSIONAL CROSSOVER IN A 
HOMOGENEOUS BOSE GAS 

In this section, we investigate the dimensional 
crossover from three to one dimension at zero temper- 
ature, as well as the dimensional crossover from three to 
two dimensions at nonzero temperature. The relevant 
geometry is illustrated in Fig. ^ 



A. Crossover from three dimensions to one 
dimension at zero temperature 

We consider a uniform Bosc gas in a box with lengths 
L x , L y , and L z and impose periodic boundary condi- 
tions. Due to the boundary conditions, the three compo- 
nents of the wave vector k take on discrete values, i.e., 
ki = 2-Kni/Li, where i = x, y, z and n, is an integer. 
Substituting this into Eq. JBJ, we obtain 



[*(*)- *(0)] a ) = E 



m 2 + \ 



i 



1 — cos 



where A^o is the number of atoms in the 
(quasi)condensate, X z = L z /2n£, and £ = 
H/[imnoT 2B (— 2/i)] 1 / 2 is the correlation length. Note 
that we have neglected the contribution from n' to the 
chemical potential, which is a good approximation in the 
weak-coupling limit we are considering. Note that the 
argument of the cosine has been simplified by choosing 
the direction of x to be along the z-axis. 

For very elongated systems, one side of the box, is 
much larger than the other two, and we choose L z 3> 
L x ,L y . If we also take L z to be much larger than the 
correlation length, we can approximate the sum over n z 
by an integral. In the limit z/^^$> 1, we obtain 



(Ix(z)-x(O)} 2 ) = ^ 



7 -_ + l n2 + ln^ 



where 



u x,y 



A, 



(9) 



(10) 



and the prime on the sum indicates that the term where 
n x = n y = is omitted. Here Kq(z) is a modified Bessel 
function of the second kind, X x y = L x y /2ir^, and 7 ~ 
0.5772 is Euler's constant. 

The one-dimensional limit is obtained by letting 
L x . y — > 0. This is equivalent to keeping only the first 
term in the right-hand side of Eq. (PJ, thus 



(IxW-x(o)] 2 ) = 



27T7Jo£ 



^+ln2 
V2 



(11) 



where uq = N$/L z is the one-dimensional quasiconden- 
sate density. The result ((111) shows that the correlation 
function g(z) at large distances falls off algebraically with 
the exponent 77 = l/47mo£- I n the weakly-interacting 
limit 47m£ 3> 1 the depletion is small, so that we indeed 
have that uq ~ n. Keeping this in mind, Eq. is 
in complete agreement with the exact result obtained by 
Haldane j2|. 

The three-dimensional limit it obtained by letting 
L XjV — > 00. The discrete sums in Eq. © then become 
integrals. Performing the integrations, we obtain in the 
limit z /£ — > 00 



V2ir 



16ir 2 



(12) 



where uq — Nq/ L x L y L z is now the three-dimensional 
quasicondensate density. Eq. i|12[l shows that the corre- 
lation function in three dimensions goes to a constant, 
so we have a "true" condensate. Note, however, that in 
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three dimensions the exponent in Eq. JSJ does not vanish 
and therefore that no appearing in Eqs. and is 
not exactly equal to the condensate density n c . In fact, 
Eq. (|12ll shows that the condensate density n c is given 
by 



n exp 



V2ir 



32tt 2 n £ 3 



8 na 3 



(13) 



where ri is the total density of the gas. This result is 
in full agreement with that of the Popov or Bogoliubov 
theory. 

The dimensional crossover behavior can be investi- 
gated in detail from Eq. © by calculating the phase 
fluctuations ~ x(0)] ) for different finite values of 

X x and X y . Comparing Eqs. (jl 1|> and l|12|l. it follows that 
the phase fluctuations scale with the product X x X y in the 
three-dimensional limit. In order to plot the phase fluc- 
tuations for different sizes of the box, it is therefore con- 
venient to multiply ([x(z) — x(0)] 2 ) by a scaling function 
f(X x ,X y ) that goes to l/X x X y in the three-dimensional 
limit. Moreover, it follows from Eq (TTJ, that f(X Xl X y ) 
must approach unity in the one-dimensinal limit. A sim- 
ple choice for a scaling function that has these properties 
is 



f(^X,Xy) — 



1 



(X X + l)(Xy + 1) 

In the following, we take X x = X y = X for simplicity. 



(14) 
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FIG. 2: Zero-temperature crossover behavior in the phase 
fluctuations of a three-dimensional homogeneous Bose- 
Einstein condensate. The three full lines show the phase fluc- 
tuations <[x(a) - x(0)] 2 )iVo/(27r) 3 (l + Xf\ z as a function of 
In (z/£) , for three values of A. The dotted lines are calculated 
from Eq. @. 

In Fig. El we plot ([ X (z) - x(0)] 2 ) Ao/ (2^) 3 (1 + X) 2 X Z 
as a function of In (z/i) for three values of A. In the three- 
dimensional limit A — > co, the phase fluctuations are in- 
dependent of z for large z/£. For finite values of A, the 



curves deviate from the ln(z/£) behavior due to contri- 
butions from the sum in Eq ©. It should be noted that 
convergence of the sum in the left-hand side of Eq. (JSJ) 
depends on the value of A. For large A, we included more 
terms in the sum to obtain convergence. The summa- 
tions were terminated when the error is less than 0.1% 
of the three-dimensional result. 



B. Crossover from three dimensions to two 
dimensions at a nonzero temperature 

The nonzero-temperature part of the phase fluctua- 
tions is given by 



- y 



N 



(15) 



The first step consists of taking the limit L y — L z — > 
oo. The sums over n y and n z then become integrals. 
The two-dimensional limit is obtained by letting L x — > 0, 
which is equivalent to keeping only the n x = term in 
the above sum. After going to polar coordinates, this 
gives 



A 2 

No Ju 



271 



d(k£)d4> 



We 



2n T 2B kiy/(ki) 2 + l/k B T 



x [1 — cos (kz cos 4>)] 
with X y — X x = X. Integration over leads to 
2ttA 2 f°° 1-J (fcz) 



Nn 



where Jo(z) is a Bessel function of the first kind. The in- 
tegral behaves as /csTln(z/^)/2noT 2B in the limit z/£^> 
1. Thus the first-order correlation function obeys 



g(z) 



l 



(z/0 1/noA2 



where A = y 2irh 2 /inksT is the thermal de Broglie wave- 
length. 

The three-dimensional limit is obtained by letting 
L x = L y = L z — > oo in Eq. i|15|) . The sums over n x , 
n y , and n z then become integrals. Going to polar coor- 
dinates and integrating over angles, we obtain 

2ttA 3 f°° d(k£) ki 1 

N Jo TP^TT ^ e 2noT^k^(k0 2 + l/k B T _ 

sin (kz) 



kz 



The first term is a constant, independent of z. The sec- 
ond term is oscillating and goes to zero in the limit 
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z/£ — > oo. This shows that the total integral, and 
thus the phase fluctuations, go to a nonzero constant 
and that Bose-Einstein condensation is possible in three 
dimensions at nonzero temperature. For low temper- 
atures, where the phonons give the main contribution 
to the above integral, this constant is proportional to 
l/no£A 2 <C 1. Combining, the above result we see that 
the crossover behavior is therefore qualitatively very sim- 
ilar to the one presented in Fig. [21 



IV. TEMPERATURE CROSSOVER IN A 
FINITE ONE-DIMENSIONAL BOSE GAS 

In this section, we first study the temperature crossover 
in a homogeneous one-dimensional Bose gas. The aim 
is to capture the essential crossover physics and to find 
the relevant parameters. Having established the rele- 
vant parameters, we examine a harmonically trapped 
one-dimensional Bose gas. 



A. Homogeneous Bose gas 

In this subsection we calculate the nonzero- 
temperature correlation function of a homogeneous 
Bose gas in one dimension. We construct the phase dia- 
gram by plotting this function versus the temperature T 
and the length L of the system. In detail the procedure 
is as follows. 

We first solve the two coupled Eqs. and J5J) for the 
(quasi)condensate density no at a fixed temperature, sys- 
tem size L, and for a total density n = N/L, where N is 
the fixed total number of atoms. Then we use this solu- 
tion in Eq. JBJ) to calculate the correlation function. To 
obtain a universal phase diagram, we need to scale T and 
L to some characteristic temperature and length in the 
system. In the present case our homogeneous system is 
characterized by its correlation length £. The correlation 
length is normally defined in terms of the condensate den- 
sity. However, above the Bose-Einstein transition tem- 
perature, the normal state cannot be characterized by 
this length. Therefore we define the correlation length 
in terms of the total density, rather than the condensate 
density, namely £ = l/\/8irK,n. As a result, we scale 
lengths to £ and temperature to To = ?i 2 /87rm£ 2 . The 
scaled temperature T/T is equivalent to (£/A) 2 , where 

A = J 2irh 2 /mk^T is the de Broglie wavelength. 

In Fig. we plot the correlation function at z — 0.9L 
versus the dimensionless variables L/£ and (£/A) 2 for 
a fixed number of atoms. This displays clearly the 
crossover from a condensate near the origin and near 
the axes, where g ~ 1, to a quasicondensate in the re- 
gions away from the origin and the axes. Fig. 0] is a 
two-dimensional contour plot that is equivalent to Fig. 13 
The contour lines in this figure show the relation between 
L/£ and (£/A) 2 for a fixed value of g. Note that when 




FIG. 3: The correlation function of a homogeneous one- 
dimensional condensate as a function of L/$ and (£/A) 2 . The 
correlation function is calculated at z = 0.9L and for a total 
number of atoms TV = 2000. 
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FIG. 4: The contour plot corresponding to Fig. The value 
of g along each contour line is indicated on the right vertical 
axis. 



quantum fluctuations can be neglected, we expect that 
-b%(0.9L)] cx (L/Ar 0(£/A) 2 (0.9L/0, which shows 
that for a fixed value of g we have that (£/A) 2 cx (L/£)~ 2 . 
This relation is indeed borne out by a detailed inspection 
of Fig. The theory discussed in this paper is only valid 
in the weakly-interacting limit Aitn£_^> 1 and therefore 
does not include the Tonks regime [25ll2c| . Nonetheless, 
it is instructive to indicate at what region of the phase 
diagram such a regime may exist. The criteria for the 
Tonks gas limit is N < AitkL [1113. Upon using the 
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above-mentioned scaling, this limit reads L/£ > y/2N. 
Since Fig. |21 was generated with N = 2000, it is clear 
that in order to be in the Tonks gas limit, the condition 
L/£> 2800 must be satisfied. 



B. Trapped Bose gas 

Next, we construct the phase diagram of a one- 
dimensional Bose gas confined by a harmonic potential, 
in order to investigate the effect of the induced density in- 
homogeneity. Our treatment is based on the local-density 
approximation for the calculation of the density profile, 
and is performed within the many-body T-matrix ap- 
proximation for the interatomic interactions (as opposed 
to the two-body limit used in the preceeding homoge- 
neous calculations). The calculation of the correlation 
function is again performed for a fixed number of atoms, 
but using now Eq. (JJJ. In the trapped case, a measure of 
the size of the system can be determined by the Thomas- 
Fermi radius Rtf, which is, in general, a function of 
both temperature and trapping frequency. To plot a uni- 
versal phase diagram, however, the axis corresponding 
to the system size should be independent of the other 
axis representing the temperature. This argument is 
strengthened by the fact that the temperature-dependent 
Thomas- Fermi radius -Rtf(T) does not increase mono- 
tonically with increasing confinement at fixed tempera- 
ture and total number of atoms. Such an effect is a direct 
consequence of two competing mechanisms which have a 
different dependence on the confinement: As the trap 
frequency increases, one would in first instance expect 
the size of the (quasi)condensate to decrease monotoni- 
cally. However, due to the constraint on the total number 
of atoms, the number of atoms in the (quasi)condensate 
tends to increase simultaneously. This can lead to an in- 
crease in Rtp(T), in particular at higher temperatures. 

To obtain a universal phase diagram in the case of 
a fixed total number of atoms, we hence use the zero- 
temperature Thomas- Fermi radius -Rtf (T = 0) as a mea- 
sure of the size of the system. This is again scaled 
to the zero-temperature correlation length £(T — 0). 
The latter is evaluated in terms o f the central density 
n(0) via £(T = 0) = 7TKn(0) = h/y/^Wjlm where 

the one-dimensional chemical potential at zero temper- 
ature in the Thomas-Fermi approximation is given by 

= (3n/V2) 2/3 (Nnl z ) 2 ^hui z |j. Note that the actual 
size of the (quasi)condensate is given by the temperature- 
dependent Rtf(T), so a good measure of the macro- 
scopic coherence of the system is obtained by measuring 
the correlation function at 0.9Rtf(T) 0- 

Fig. shows the value of the normalized correlation 
function of Eq. JSJ evaluated at z — 0.9Rtf(T) as a 
function oiR TF (T = 0)/£(T = 0) and (£(T = 0)/A) 2 for 
fixed coupling constant k and total number of atoms N in 
the system. Light color indicates the presence of a pure 
condensate and black color indicates the absence of uni- 




FIG. 5: Universal phase diagram of a pure one-dimensional 
trapped atomic Bose gas with constant number of atoms. 
Greyscale indicates the value of g(0,z — 0.9Rtf(T)) with 
lighter regions corresponding to maximum coherence g ~ 1. 
This figure probes the entire region from the normal-to- 
(quasi)condensate crossover (where the (quasi)condensate de- 
pletion approaches 100%, curved region away from origin), to 
the presence of true Bose-Einstein condensation, for which 
the quasi-condensate depletion becomes negligible near the 
origin. The entire phase diagram plotted here remains well 
within the validity of our Thomas-Fermi approximation, with 
7 = NqkIz S> 1 and iVo the number of (quasi)condensate 
atoms. For a 23 N a condensate with N = 2000 atoms under 
very tight transverse confinement of uj±/2tt = 5 KHz, the 
above phase diagram would correspond to temperatures in 
the range 0.2nK — 12fiK and longitudinal frequencies ui z /2n 
from 0.05 to 45000 Hz. 



versal phase coherence across the system. As expected, 
the coherence is found to increase when decreasing the 
size of the system, i.e., increasing the trap frequency uj z 
at fixed temperature, or decreasing the temperature in 
a fixed trap. This is similar to the homogeneous case, 
with g becoming close to one as Rtf(T = 0)/£(T = 0) 
or (£(T = 0)/A) 2 become small. Note that the quasicon- 
densate depletion due to interactions and nonzero tem- 
peratures can be very significant, and this figure covers 
the entire range of (quasi)condensate fractions, ranging 
from one at the origin, where we essentially deal with a 
pure condensate, to approximately zero near the transi- 
tion to the normal state sufficiently far from the origin. In 
particular, the characteristic curved shape away from the 
origin marks a sharp transition from a quasi-condensate 
to the normal phase. Note that we find a sharp transi- 
tion due to the use of the local-density approximation, 
whereas in reality this would be a smooth crossover. To 
study now the crossover between regions of different de- 
gree of coherence within the (quasi)condensate regime, 
Fig. EJa) shows the corresponding contour plot, focusing 
on a small region near the origin, where the crossover 
takes place. In this figure, white regions correspond to 
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true Bose-Einstein condensation, whereas the darkest re- 
gion on the top right corner of the figure corresponds to 
complete absence of phase coherence across the size of 
the system. The corresponding contour lines with both 
axes plotted on logarithmic scale is shown in Fig.EJb). 

(a) 



0.25 




50 100 t50 



R TF (T=0)/^(T=0) 
(b) 
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FIG. 6: Contour plot corresponding to Fig. O with both axes 
plotted on (a) linear and (b) logarithmic scale. As in Fig. 2] 
each contour corresponds to a value of g ranging from 0.1 
(leftmost contour) to 0.9 (rightmost contour). White regions 
correspond to true Bose-Einstein condensation, and darkest 
regions to the absence of phase coherence over the system size. 
Note that greyscale in this figure is not the same as in Fig. 
Note also that we have not explicitly computed for Rtf(T = 
0)/£(T = 0) < 7, where there is essentially true Bose-Einstein 
condensation, as this requires extremely high temperatures 
and large longitudinal trapping frequencies. The condition of 
validity of the Thomas-Fermi approximation will also break 
down in this limit. 

The above phase diagram corresponds to the limit of 
a pure ID Bose gas. In trying to probe this experimen- 
tally by means of the recently produced quasi-lD Bose 



gases 0, 0, 0, 0, there are certain constraints 
which must be taken into consideration. First, the tem- 
perature must be low enough that transverse thermal 
excitations are suppressed, i.e., kT < 7kJj_. Second, the 
chemical potential must not exceed the transverse energy, 
i.e., /i < hu)± thus restricting the maximum longitudi- 
nal frequencies for which the system remains kinemati- 
cally one dimensional, to some fraction of the transverse 
confinement. The exact fraction depends on the atomic 
species, the scattering length, the number of atoms and 
the temperature. Finally, the longitudinal confinement 
cannot be made arbitrarily weak, as it will then be es- 
sentially impossible to trap the atoms. All above con- 
straints require a very large transverse confinement, since 
the experimentally-accessible region of such a phase dia- 
gram increases with increasing transverse confinement. 
The experimental constraint of minimum longitudinal 
frequency does not pose much of a limitation since the 
low frequency regions correspond to the outmost points 
of Fig. |S1 for which the macroscopic coherence has already 
been completely lost. Of the other two constraints, which 
is more restrictive depends on the details of the partic- 
ular experiment, since they exclude different regions of 
the phase diagram. 

To illustrate this, Fig. 0[a) indicates the experimen- 
tally accessible regions for the case of N = 20 23 Na 
condensate atoms under a transverse confinement of 
uj±/2-k = 5 KHz j2^. The chemical potential constraint 
generally excludes a region to the left of Fig.|SJa), and so 
large coherences are best probed by restricting the con- 
densate to a very small number of atoms. For the case 
studied here, this translates to the approximate condition 
R TF (T = 0)/£(T = 0) > 40. The temperature constraint 
on the other hand prohibits a large part of the top left of 
Fig-E^a), which requires the experiments to be performed 
at sufficiently low temperatures (here T « 2A{)nK). We 
see that for low enough temperatures and small enough 
atom numbers, it is possible to probe experimentally the 
entire phase coherence regime (from g = to roughly 
g = 0.99). Although the coherence decreases with in- 
creasing scattering length, for example by tuning around 
a Feshbach resonance [23, H3]> or central density, this 
can further restrict the region of accessibility due to the 
chemical potential constraint. 

Fig. 7(b) gives the corresponding contour plot of the 
(quasi)condensate fraction. It is evident that the 
(quasi)condensate depletion becomes important even 
within the range of this figure, with the furthermost 
points here corresponding to more than 50% depletion. 
In our analysis, we have chosen to monitor the loss of 
phase coherence by looking at the amount of decrease of 
the normalized two-point correlation function near the 
edge of the system (at z = 0.9Rtf )■ However, it may 
be experimentally easier to look at the decay of this cor- 
relation function at a point somewhat closer to the cen- 
tre. This would then lead to a slower decay of phase co- 
herence across the measured distance, and a sufficiently 
small amount of phase coherence across such a region 
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FIG. 7: (a) Same contour plot as in Fig. [(J a) indicating clearly 
the regions experimentally accessible for N — 2, 00 23 Na 
atoms with u>±/2tt — 5 KHz. The solid line corresponds to 
the condition fj, — hu± , whereas the dashed line indicates the 
regime where kT = hui±. Experimentally accessible regions 
due to actual quasi-lD constraints are located to the right of 
each of these lines, (b) Corresponding contour plot of quasi- 
condensate fraction. Each contour corresponds to a change in 
No/N = 0.05, with the leftmost contour indicating the region 
where No/N = 0.95, and the rightmost contour in the plot- 
ted region corresponding to No/N = 0.5. It is thus evident 
that the quasi-condensate depletion becomes significant very 
rapidly, as one moves away from the origin of the figure. 



would then occur for even larger (quasi)condensate de- 



pletions. Therefore, we believe that a quantitative study 
of the phase-coherence crossover problem requires a the- 
ory which takes into account this depletion. 

V. CONCLUSIONS 

We have obtained a detailed universal picture of the 
coherence properties of a Bose g function of 

the dimensionality, the temperature and the interaction 
strength, paying particular attention to the crossover 
from quasicondensation to true Bose-Einstein condensa- 
tion. In a finite homogeneous system at zero tempera- 
ture, the long-range coherence is lost when the dimension 
of the system is reduced to one dimension. This dimen- 
sional reduction leads to the appearance of a quasicon- 
densate, whose coherence extends over a range that is 
much larger than the coherence length but much smaller 
than the size of the system. Similar results were ob- 
tained when the effects of temperature and the strength 
of the mean-field interaction were investigated. In a ho- 
mogeneous one-dimensional system, a true condensate is 
located in a region in the phase diagram where both the 
length of the system and its temperature are small. Away 
from this region, a quasicondensate is present covering a 
wider region of the phase diagram. 

A similar phase diagram was obtained for a harmoni- 
cally confined system in one dimension, where a change 
in the confining potential at fixed temperature leads to a 
change in the size of the (quasi)condensate. It is impor- 
tant to note that our theory covers the entire temperature 
range, from zero temperature up to the crossover to the 
normal phase, where the (quasi)condcnsate becomes fully 
depleted. Since our work takes the (quasi)condensate 
depletion into account, the presented results extend be- 
yond existing treatments which have been limited to the 
regime of negligible depletion 0, 0, ^| . Note also that 
the phase diagram presented in Sec. IV corresponds to 
the pure one-dimensional limit. Nonetheless, a large part 
of this, ranging from, complete absence of, to complete 
coherence over the system, can be probed experimentally 
with existing technology, as discussed by a particular ex- 
ample based on typical parameters. 
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